Poisson's ratio in composite elastic media with rigid rods 



O 

(N 

< 



Moumita Das and F.C. MacKintosh 
Department of Physics and Astronomy, 
Vrije Universiteit, Amsterdam, The Netherlands. 

(Dated: September 1, 2010) 

We calculate both the micro-mechanical response and bulk elastic constants of composites of rods 
embedded in elastic media. We find two fixed points for Poisson's ratio with respect to rod density: 
there is an unstable fixed point for Poisson's ratio=l/2 (an incompressible system) and a stable fixed 
point for Poisson's ratio=l/4 (a compressible system). We also derive approximate expressions for 
the elastic constants for arbitrary rod density, which agree with exact results for both low and high 
density. These results may help to explain recent experiments [Physical Review Letters 102, 188303 
(2009)] that reported compressibility for composites of microtubules in F-actin networks. 
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There are many natural examples of composite materi- 
als, combining multiple components with distinct elastic 
properties. Wood, bone and various tissues are all made 
of composites, as is each living cell [H. Many composites 
consist of a reinforcing constituent such as stiff fibers em- 
bedded in a weaker, less stiff matrix Q. By varying the 
relative concentration of the constituents, one can tune 
the elastic moduli, resulting in materials with remarkable 
properties 0, Q ■ The collective properties are more than 
merely the sum total of those of its constituents. Among 
the most prevalent natural or living composites is the cell 
cytoskeleton, which consists of a complex scaffold of sev- 
eral distinct filamentous proteins, some of which are very 
rigid. Most previous biophysical studies of cytoskeletal 
networks have focused on purified gels or networks con- 
sisting of one type of filament 5-14]. The cytoskeleton, 
however, contains three major types of filaments: micro- 
tubules (MTs), filamentous actin (F-actin), and inter- 
mediate filaments. These filaments have vastly different 
bending stiffness. A few studies of reconstituted compos- 
ite cytoskeletal networks have shown viscoelastic proper- 
ties distinct from single-component networks jlSl - llTj . A 
very recent experimental report also provides evidence 
of anomalous compressibility with the addition of stiff 
microtubules to a soft matrix [Tsj . 

Here, we develop a model for the mechanical response 
of a composite material consisting of rods in an elastic 
matrix, using a mean field approach and a dipole approx- 
imation for the rod-like inclusions. The elastic matrix 
under consideration is treated as an effective medium 
that is made of the bare elastic medium (e.g., the F- 
actin matrix) and a collection of rods (MTs) embedded in 
it. Consistent with the experiments of Ref. [l^, we find 
that the addition of rigid rods can lead to enhanced com- 
pressibility of an initially nearly incompressible medium. 
Specifically, we find that for matrices characterized by 
Poisson's ratio 1/4 < i/ < 1/2, the addition of rods re- 
duces I', while for i' < 1/4, stiff rods increase i'. In this 
way, — 1/4 can be thought of as a stable fixed-point 
of such a composite. We further evaluate the Poisson's 
ratio and elastic moduli as functions of the concentration 



of rod-like inclusions. While this mean-field approach is 
only approximate at intermediate concentrations, we ob- 
tain an exact result in the limit of high concentration. 

We first study the micro-mechanics of our system using 
the elastic response function or Green's function for the 
displacement field u in response to an applied force. For 
an isotropic and homogenous elastic material with Lame 
coefficients A and /i, we can describe the displacement 
field Ui at a position r in the medium due to a force / 
acting at point r', Ui(r) — aij{f~r')fj{r'), with the linear 
response function aij given by : 
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where /3 = ^/{X + 2/z) is the ratio of the shear modu- 
lus fi to the longitudinal modulus. For an isotropic and 
homogeneous elastic material, Q:ij{f) reduces to just two 
distinct components corresponding to the response par- 
allel and perpendicular to r, as shown in Fig. [T^. For 
an incompressible material in 3d, P — and the paral- 
lel and perpendicular response functions are related by a 
simple factor of two: Q:\\{r) = 2a_i{r) = 1 / 47r u r, which 
is the elastic analogue of the Oseen Tensor [l^ [20| • For 
compressible elastic media in 3d, /3 > 0, a||(r) = l/Anfir 
and a±{r) = (1 -|- 13) /Sir fir. Thus, /3 provides a measure 
of the degree of compressibility. 

We calculate the change in the response function and 
Lame coefficients upon addition of rods as follows. We 
consider a single rod of length a embedded in the elastic 
medium, as shown in Fig. [1] (a). The presence of the rod 
represents a constraint on the displacement field induced 
by the applied force. For a force fy applied at A, the net 
displacement of the ends of the rod of length a oriented as 
in Fig.[T](a) is given by AM(r) = {u{f+a/2) — u{f—a/2)). 
We approximate the constraint of an incompressible rod 
by a dipole at its center of mass. This induced (tensile) 
dipole is oriented along the rod and its strength is chosen 
so as to enforce a constant end-to-end distance of the 
rod: p — fnra{d ■ Au)/2 — x^, where x = I^T^a^ /2, e = 
a ■ (a ■ V)u, and where we keep only leading-order terms 
in a, which is assumed to be small compared to the other 
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FIG. 1: (Color Online) We consider a point-force jy applied 
at the origin (^) and calculate the response at points B lo- 
cated parallel (here, along the y— axis) and perpendicular 
(along the 2— axis) to the applied force in Fig. 1 (a). The 
rod center of mass is located at a distance r from the origin, 
and makes a polar angle Q and azimuthal angle 0. The rod is 
oriented in a direction as shown in Fig. 1 (b). 



lengths in the figure. The resulting displacement field at 
B allows us to determine the change in the linear response 
functions. In our effective medium picture, we obtain 
the change in response as arising from a cloud of induced 
dipoles in the elastic continuum. 

The change in the parallel and perpendicular response 
functions with the addition of the rods are 6a\\ — 
— (7r/30)na^a|| and (5ckj^ = ^(1 -I- 3/3^)5aj[, where n is 
the rod number density and we have averaged over rod 
orientation. This is equivalent to 5\i = 8\ = j^Xf^- For & 
small increment dn in added rods, we obtain the differen- 
tial equation d(3/dn — j^a'^/3(l — 3/3). Thus, /3 increases 
for < /3 < 1/3 and decreases for (3 > 1/3, while f3 
is unchanged for /3 = and (3 — 1/3. Therefore, for a 
slightly compressible medium, adding rods enhances the 
compressibility relative to the shear compliance, while 
for a highly compressible medium, the rods reduce the 
compressibility, as illustrated in Fig. [21 For finite com- 
pressibility, addition of rods tends to drive the system 
towards a state with (3 = 1/3. This suggests a stable 
fixed point (to the addition of rods) at (3 — 1/3 and 
hence Poisson's ratio v — 1/4. Similarly, /3 = and 
1/ — 1/2 corresponds to an unstable fixed point. These 
results are also apparent from the full solution to the dif- 
ferential equation above for /3, which is shown in Fig. [3] 
for several different initial values /3o in the absence of 
added rods. Our results are qualitatively consistent with 
recent microrheology experiments on a coinposite of mi- 
crotubules embedded in filamentous actin [l^ , which re- 
ported enhanced compressibility {v < 0.5) when stiff mi- 
crotubules were added to an almost incompressible actin 
matrix {ly ~ 0.5), as inferred from the measured parallel 
and perpendicular response functions. 

We can also calculate the effect of the addition of stiff 
rods by considering uniform strain of a composite, as 
follows. We apply a uniaxial strain ezz to the medium 
along the z direction and clamp its boundaries along the 
X and y directions. We consider a particular rod making 
polar and azimuthal angles 6 and with respect to the 
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FIG. 2: (Color Online) The flow diagram for the degree of 
compressibility (3 showing a stable fixed point at /? = 1/3, 
and an unstable fixed point at /3 = 0. 



coordinate shown in the schematic Fig. [IJb). In 

a way similar to the above, the presence of this rigid 
rod is equivalent to a dipole, here of strength p — xe, 
where x is the same as above, and e — t^z cos^ d. We 
calculate the additional stresses on the boundaries due 
to an isotropic distribution of such rods with number 
density n. These stresses are given by ^a^x = SXczz and 
Sczz = {26^ + 5X)ezz- For rod orientations in a given 
solid angle dil = sm 9d9d(p, the stress associated with 
the induced dipole is given by Satj = j^x^^iO-jd^- Thus, 
the Lame constants can be calculated by 

SX ^ nx f cos(6l)2sin(6')2cos(0)2 — 
J 47r 



{2Sn + SX) = nx / cos(6l) 
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Solving for 6fj, and i5A, we find the same values as in the 
previous micromechanical calculation. 

Wc now consider the rods to have finite stretch mod- 
ulus K ~ nb'^Er, where b is the radius of the rod with 
Young's modulus Er- When the rod is subject to an ex- 
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FIG. 3: The degree of compressibility f5 and Poisson's ratio u 
of the composite as a function of mesh size ^ for inextensible 
rods, for different values of the degree of compressibility /3o of 
the medium in the absence of rods. The mesh size ^ is related 
to the rod density nhy 1/^^ = na. 
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tensional strain e, the resulting force balance for an ex- 
tension A of the rod is now eiven by A = ae = ^ 

where p = xe, and x — TTa^^/{2 + ira^^/K). For a den- 
sity n of rods, we obtain the following nonlinear relations 
for the Lame coefficients: 



S\ = 6fi = ^^^na^/ ( 1 
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In the limit of very stiff rods, this reduces to the expres- 
sion above. For highly compliant rods, by contrast, this 
is consistent with the elastic moduli of an affinely de- 
forming rod network of volume fraction (p — nab^n; e.g.. 

The result in Eq. ([3]) is valid at small density, in which 
the shear modulus fi on the right-hand-side is that of the 
(bare) matrix. If we consider n = dn to be a small in- 
crease in the number density of rods, then Eq. ^ can 
be thought of as a set of differential equations repre- 
senting the increase of the moduli upon the addition of 
stiff rods. By treating the resulting composite system 
as an isotropic and homogeneous effective medium, this 
differential equation suggests a way of calculating the 
properties of composites with finite rod density. This is 
similar to self-consistent methods employed for aligned 
fiber-reinforced composites [21j|. While this represents 
an uncontrolled approximation, we find that integrating 
Eq. ^ yields an exact expression for the limit of a high 
density of rods. The solution for fi is given by 



H ^ /.irW — exp 
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where, fiQ is the shear modulus of the medium in the ab- 
sence of rods and /i,- = 2K/ (tto^). Here, W(z) is the prin- 
cipal value of the Lambert W-function, which is defined 
hy z = We^ . From this, we also obtain X = Xo + fi — fio- 
The results are shown in Fig. |3] for various initial condi- 
tions fj,o and Aq. For small n and large K, this reduces 



to ^ ~ ^o(l + 7rna''/30), since W{z) 



for small 



This is consistent with the above results for dilute, in- 
extensible rods. As the density of rods and correspond- 
ing shear modulus increase, however, W{z) ~ ln(z) and 
/i — > j^(pEr, which is the modulus of a high density mesh- 
work of elastic rods. 

So far, we have not accounted for the tension profile 
along the rod. Our dipole approximation is expected to 
overestimate the effect of the rod, since the real displace- 
ment field along a finite rod is expected to vary more 
smoothly than for a dipole. The displacement field cor- 
responds to a strain and consequent tension along the 
rod that is uniform at its center and vanishes at its ends. 
This strain field can be calculated from the force balance: 
Kv"{x) — C, [v(x) — ex], where v(x) is the displacement 
field along the rod in the presence of a background strain 
e of the matrix, and C, represents the elastic coupling of 
the rod to the matrix. The longitudinal strain of the rod 
is given by v'{x), and the gradient of this corresponds to 
a net force per unit length on the rod, which we take to 
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FIG. 4; (Color Online) The solid blue lines show the shear 
modulus of the composite as a function of the mesh size for 
different values of the ratio of the rod and medium compli- 
ance, while the dashed red line represents the afiine result. 
The inset shows the Lame coefficient A for different values of 
the initial degree of compressibility /3o of the medium for both 
extensible and inextensible rods. The ratio fio/fir ~ 0.001 in 
the inset, except for the inextensible rod (dashed line), where 
/3o = 0.1. 



be proportional to the displacement of that section of the 
rod relative to the background medium. This is similar 
to the viscous drag on a slender body in the presence of 
a background velocity field (2^], in which case C is the 
drag coefficient per unit length of the rod. For an elastic 
medium, we approximate ( = 27r/i/ ln(^/6), as in Ref. 
[23 |. Here, the screening length ^ is taken to be of order 
the average separation or mesh size of the rod network, 
which depends on density. Given the weak logarithmic 
dependence on ^, however, we will treat ^ as a constant. 

Once again, the rod can be considered as a force dipole 
on scales large compared with the rod length a. Using 
the condition of tension free ends, the strength of this 
dipole can again be expressed as p = but now with 



X = Ka[l- 24 tanh {a/2£o) /a] 



(5) 



Here, Iq = yK/C represents the length over which the 
longitudinal state of strain of the rod varies J^S]. For 
highly compliant rods, this becomes a small length, cor- 
responding to a nearly constant state of strain and ten- 
sion along the rod, except very close to the ends. In 
the other limit, of very stiff rods, the strain exhibits a 
quadratic dependence, with a maximum at the center 
of the rod and vanishing at the ends of the rod. In 
this case, x = C^^^/l^, which is smaller by a factor of 
31n(^/&) than the value above for the simple dipole ap- 
proximation. Interestingly, independent of the parame- 
ters of the system, we still find that both Lame coeffi- 
cients evolve in the same way upon the addition of rods: 
dfj, = dX = ?T.x/15 [i^l. This means that the qualitative 
form of dp/dn in Fig. [21 as well as our conclusions re- 
garding the fixed points at /3 = and j3 = 1/3, remain 
unchanged. For rods that only interact with each other 
through their matrix (either a viscous fluid or a homo- 
geneous elastic matrix), the concentration of rods only 



4 



enters the calculation via the screening length ^. Fur- 
thermore, the modulus appearing in the coefficient ^ is 
that of the bare matrix. For rods that interact directly 
with each other, we employ a self-consistent approxima- 
tion, as discussed above [2l| . A more detailed analysis of 
the effects of direct inter-fiber interactions would likely 
require a numerical simulation that is beyond the scope 
of the present work. Nevertheless, we can derive a dif- 
ferential equation for /i accounting for the tension profile 
along the rods using our self-consistent approach. Al- 
though this is more complicated than Eq. ^ , the dipole 
strength derived from Eq. ([5]) is well approximated by 
that used in Eq. ([3]), apart from the a factor of 3 In (^/b). 
Thus, the functional forms in Figs. [3] and |4] are expected 
to be good approximations. 

We have studied the collective mechanical response of 
composites of rods embedded in elastic media, such as 
stiff MTs in a softer cytoskeletal matrix or carbon nan- 
otubes in synthetic gels ^25J , using a mean field approach 
similar in concept to homogenization methods introduced 
for elastic composites . We find a very general result 
that the addition of elastic rods or fibers leads to a mono- 



tonic evolution of Poisson's ratio toward the value 1/4, 
either from above or below. On the one hand, this is 
consistent with recent numerical calculations for fiber- 
reinforced concrete, showing a weak increase in v with 
fiber density in the range 0.2 < v < 0.25 [1^. On the 
other hand, our results may help to explain recent exper- 
iments [18] that have reported < 1/2 for composites 
of microtubules and F-actin networks, while u ^ 1/2 
for single-component F-actin networks. This suggest an 
important role for stiff filaments such as MTs and stress 
fibers in the mechanics of the cell cytoskeleton — they not 
only enhance the stiffness of the cytoskeleton [l^] and 
its ability to bear large forces, but may also endow it 
with enhanced compressibility relative to the shear com- 
pliance. 
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